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SECOND-ORDER SUBSONIC AIRFOIL THEORY INCLUDING EDGE EFFECTS ^ 

By Milton D. Van Dtxb 


SUMMARY 

Several recent advances in 'plane subsonic jlow theory are 
combined into a unified second-order theory for airfoil sections 
of arbitrary shape. The solution is reached in three steps: 
The incompressible result is found by integration, it is corv- 
verted into the corresponding subsonic compressible result by 
means of the second-order compressibility rule, and it is ren- 
dered uniformly valid 'near stagnation points by f'urther rules. 
Solutions for a number of airfoils are given and are compared 
with the results of other theories and of experiment. A straight- 
forward computing scheme is outlined for calculating the sur- 
face velocities and pressures on any airfoil at any angle of 
attack. 

INTRODUCTION 

Thin-airfoil theory provides a useful first approximation 
to the incompressible flow past two-dimensional airfoils, 
and the results can be immediately extended to subsonic 
compressible flow by the Prandtl-Glauert rule. It is natural 
to attempt to improve this simple theory by successive 
approximations so as to increase its accuracy for thicker 
airfoils and higher subsonic Mach numbers. There results 
a series expansion of the flow quantities in powers (supple- 
mented by logarithms in the fourth and higher approxima- 
tions) of the airfoil thickness ratio, camber ratio, and angle 
of attack. 

For incompressible flow, the higher-order theory has been 
studied by various writers, in particular Riegels and Wittich 
(refs. 1 and 2) and Keime (ref. 3). A less straightforward 
series of approximations was developed by Goldstein (ref. 
4). Perhaps the most concise exposition of higher-order 
incompressible thin-airfoil theory is given by Lighthill 
(ref. 6). 

For subsonic compressible flow, the corresponding analysis 
was first undertaken by Gbrtler (ref. 6), followed by Hantzsche 
and Wendt (refs. 7 and 8), Schmieden and KawaDd (ref. 9), 
Kaplan (refs. 10 and 11), Imai and Oyama (refs. 12 and 13) 
and others. These investigators treated only specific simple 
shapes by rather laborious analysis. Later, it was discovered 
that particular integrals of the second-order iteration equa- 
tion can be expressed in terms of the first approximation 
(refs. 14 and 15). This permits the second-order subsonic 
solution for any profile to be given in terms of integrals 
(refs. 16 and 16). However, the resulting solutions are 
sometimes incorrect everywhere for airfoils with stagnation 
points, for reasons to be discussed later. 


Recently Hayes (ref. 17), improving on a result of Imai 
(ref. 18), has given a second-order similarity rule for surface 
pressure that implies a second-order extension of the Prandtl- 
Glauert rule (ref. 19) . This remarkable result was overlpoked 
by earher investigators because they did not calculate sur- 
face pressures, but were content with finding surface speeds, 
for which the second-order compressibility rule is more 
complicated. These rules reduce the second-order problem 
of subsonic compressible flow past airfoils to the correspond- 
ing incompressible problem. 

However, the solution by successive approximations 
breaks down near leading and trailing edges if there are 
stagnation points. The result is therefore merely a formal 
series expansion, which fails to converge near the edges. 
In first-order theory spurious singularities arise at stagnation 
edges, but it is known how they can be taken into account, 
since they are integrable. In the second approximation, 
however, these singularities are intensified, so that at round 
edges they are no longer integrable. In any case, the calcu- 
lated speeds and pressmres are incorrect near such edges, 
and more so in the second approximation than the first. 
Moreover, in subsonic compressible flow the second ap- 
proximation may be incorrect everywhere as a consequence 
of the defects in the first approximation. 

For roimd edges in incompressible flow, previous in- 
vestigators have shown hoV these defects can be corrected. 
Riegels (ref. 2) gave a simple rule that renders the first- 
order thin-airfoil solution valid near the edge. LighthiU 
(ref. 5) gave an equivalent rule for the second approxima- 
tion. Recently, corresponding rules have been developed 
for higher approximations, for sharp as well as round edges, 
and for subsonic compressible flow (ref. 20). 

It is the aim of this paper to combine these recent advances 
into a unified theory. There results a unif orm second ap- 
proximation to subsonic flow past any profile at angle of 
attack, expressed in terms of integrals that can, if necessary, 
be evaluated nmnericaUy. It may be noted that, except 
possibly for certain particular shapes at isolated Mach 
numbers, the resulting solution is now generally believed 
to be valid only below the critical Mach munber — that is, 
for pmely subsonic flows. Although only flow quantities 
at the airfoil surface are considered here in detail, the entire 
flow field could be treated in the same way. 

For numerical computation, the most useful method 
appears to be that initiated by Riegels and Wittich (refs. 
1 and 2) and independently by Germain (ref. 21), with 


> Supersedes NAOA TN 3390, “Second-Order Subsonic AlrfoU-Sectlon Theory and Its Practical Application,” 1955, and portions ol NAOA TN 3343, “Subsonic Edges In Thln-Wlng and 
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extensions by Watson (ref. 22), Thwaites (ref. 23), and 
Weber (ref. 24 and 25). It requires a knowledge only of 
the airfoil ordinates at a specified set of points. In this 
report a straightforward scheme, based on an extension of 
this method, is given for computing the second-order sub- 
sonic solution for any airfoil. The reader interested only 
in calcidating a specific case, without necessarily under- 
standing the details of the theory, can turn directly to the 
section “PRACTICAL NUMERICAL COMPUTATION.” 

The author is indebted to R. T. Jones for many helpful 
discussions throughout the course of this work. 

THEORY 

From the preceding remarks it is clear that the solution 
is reached in three steps. First, the formal second-order 
incompressible solution is foimd by integration. Second,, 
this is converted into the corresponding subsonic com- 
pressible solution by means of the second-order compressi- 
bihty rule. Third, this is modified near stagnation points 
by the appropriate rules for roimd or sharp edges. These 
three steps wiU be considered successively. 

FORMAL INCOMPRESSIBLE SOLUTION 

The expansion of the velocity components in a formal 
series of powers of the airfoil thickness ratio, camber ratio, 
and an^e of attack has been discussed in detail by Lighthill 
(ref. 5). It win suffice here to summarize his results for the 
second approximation. We mainly foUow his notation 
except to make it more mnemonic, and to suppress his 
parameter e characteristic of the airfoil thickness, which is 
only convenient in the detailed analysis. 

Accordingly, consider an airfod of moderate thickness and 
camber at a moderate angle of attack to a uniform sub- 
sonic stream (fig. 1). It is essential that the x axis be 
chosen to pass through both the leading and trailing edges. 
Let the upper and lower surfaces of the airfoil be described by 

2/=F(a;) = (7(a;)±r(a:) (1) 


y 



The corresponding horizontal velocity disturbance on the 
chord line, which is required for calcidating the surface 
pressure, consists of a term associated with the airfoil 
thickness, and another associated ivith its camber and 
angle of attack. For the thickness 

^=1 (f (3) 

U irjia x—i 


and for the camber and angle of attack 
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The latter result is due to Munk (ref. 26) and tlie former 
was first given by Pistolesi (ref. 27). Cauchy principal 
values are indicated in each integral. 

The surface speed is then given to a first approximation by 


Qi I t I '^ic 


( 6 ) 


Second-order solution. — In the second approximation, the 
tangency condition is transferred from the airfoil surface 
to the chord line by Taylor series expansion. The condition 
on the second-order increment in vertical velocity is thus 
found to be 


where 



C',(x)±T',(x) 


a(x)=^ T 

T,(x)=^ O 


(6a) 


(6b) 


(We depart here from Lighthdl’s notation in order to em- 
phasize that the functions Ci and 2j are effectively the 
camber and thickness for some fictitious airfoil.) The 
problem is identical with that in first-order theory except 
for the condition at infini ty, which is readily disposed of. 
Thus, corresponding to T 2 is the increment in horizontal 
velocity 


«?£_1 1^ 
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and corresponding to C 2 
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Figure 1. — Notation for airfoiL 

where C{x) describes the mean camber line and T(r) the 
thickness. The airfoil extends over the interval Xa<x<Xj,, 
which is usually conveniently taken to be either — 1 <a:<l 
or 0<a:<l. All symbols are defined in Appendix A. 

First-order solution. — In the first approximation of thin- 
airfoil theory, the condition of tangent flow at the airfoil 
surface is imposed on the two sides of the chord line y=0 
rather than at the surface, and requires that 


The velocity components on the surface of the airfoil 
include also terms arising from the transfer from the chord 
line to the surface, which is again effected by Taylor series 
expansion. Hence the surface speed is given to a second 
approximation by 

g=H-^±|)?+^±^+(<7±D(C7”±T”)-fi {0'±T'y (0) 
and the surface pressure coefficient by 
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Higher-order solutions can be found by continuing this 
process. Lighthill gives explicit formulas for the third-order 
solution. 

Airfoil integrals. — ^The incompressible solution to second 
order (or, indeed, to any order) is thus reduced to a succession 
of “airfoil integrals” typified by equations (3), (4), (7), and 
(8) . Goldstein (ref. 28) emphasizes that in first-order theory 
these integrals can be evaluated analytically for practically 
every profile for which formulas have ever been proposed. 
In second-order theory this appears to true to a somewhat 
lesser extent, although the labor of calculation becomes great 
except for simple shapes. Often the integrals are most 
readily evaluated by guessing (u—iv) as a function of the 
complex variable {x-\-iy) that has the required behavior on 
the chord line. A short table of airfoil integrals useful for 
finding second-order solutions is given in Appendix: B. 
Other can be found in references 28, 29, and 30. 

For complicated profiles, exact analytic evaluation of the 
integrals may be impossible or excessively laborious. Then 
munerical integration may be resorted to, or the profile can 
bo approximated by a simpler shape that can he treated 
analytically. The most useful numerical procedme is ap- 
parently that originated independently by Riegels and 
Germain and simplified and extended by Watson, Thwaites, 
and Weber. In this method the airfoil ordinates are approxi- 
mated by the trigonometric polynomial 

N-l 

Y (^r cos rO+tr sin r6)+k{f cos Nd 

r=l 




Fiqueb 2 . — Parametric angle d. 



whore 9 is the angle indicated in figme 2. The coeflficients 
kr (for camber) and tr (for thickness) are chosen to give the 
actual ordinates at the 2N points for which 6=jmrlN. In 
this way it is found that the airfoil integrals can be expressed 
approximately as s ums of the airfoil ordinates at certain 
pivotal points multiplied by standard influence coefficients. 
The details of this method, as adapted to second-order thin-air- 
foU theory, are given in Appendix C. The numerical com- 
puting procedure is outlined in the last section of this paper. 


SECOND-ORDER COMPRESaiBILITT RULE 

The second-order coimterpart of the Prandtl-Glauert 
compressibility rule is implicit in an extension of transonic 
s imili tude that was initiated by Imai (ref. 18) and carried 
to completion by Hayes (ref. 17). Imai sought to improve 
the transonic s imil arity rule by retaining in its derivation 
aU terms proportional to the square of the airfoil thickness 
except one appearing in the condition of tangent flow at the 
surface. The correlation of experimental data was not 
appreciably unproved, which led him to surest that the 
neglected second-power term should also be included. This 
probably cannot be done for the whole flow field. However, 
in attempting merely to reproduce Tmai’s result as announced 
before publication, Hayes actually included that term in a 
second-order rule for surface pressure. 

Hayes’ result is that for two-dimensional subsonic or super- 
sonic flow the ratio of the second-order to first-order pressme 
term on the surface is proportional to the parameter 

gj3^[5±lM*+2a-ilO] (ID 

where t is some measure of the thickness, camber, or angle 
of attack. Now at subsonic speeds the first-order pressure 
term is related to its value in incompressible flow by the 
Prandtl-Glauert rule. Combining these two results yields 
the second-order compressibility rule (ref. 19). 

In incompressible flow the second-order surface-pressme 
coefficient has the form 


C^ix) = C^ix)+AC^ix) (12a) 

where the first-order term Cp^ contains linear terms in thick- 
ness, (jamber, and angle of attack, and the second-order 
increment AOp^ contains their squares and products. Then 
for the same airfoil in subsonic compressible flow, according 
to the compressibility rule, the pressure coefficient is 


where 


1 1 


Ki= 




{y+l)M*+40^ 


(12b) 


(12c) 


The corresponding compressibility rule for surface speed 
is readily foimd from the above rule for pressure by con- 
sidering the smaU-disturbance series form of BemoiiUi’s 
equation for compressible flow. Thus it is found that if 
the surface speed ratio in incompressible flow is 


So. 
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where Agi contains linear terms in thickness, camber, and 
angle of attack, and A23 their squares and products, then 
at subsonic speeds 
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2 ' 8j3< 


(13c) 


This rule is seen, to lack the fundamental simphcity of the 
rule for pressme. 


EDGE COEBECnONS 

Thin-airfoil theory is known to fail near leading and 
trading edges if there is a stagnation point. The flow is 
actually brought to rest, hut thin-airfoil theory predicts . 
infinite speeds instead. If r is the distance from the edge, 
the velocity contains powers of r-1/2 for a roimd edge and 
for any leading edge with flow arovmd it (associated with 
angle of attack), and powers of Inr for a shaip edge. First- 
order theory contains first powers of these singularities, 
second-order theory their squares and products, and so on, 
so that the formal thin-airfoil series diverges in some neigh- 
borhood of the edge. Not only are the velocities and pressure 
incorrect near stagnation edges, but nonintegrable singu- 
larities appear in the higher-order e.xpressions for aero- 
dyanamic forces on round edges. 

False subsonic solutions. — ^Even more serious difficulties 
may arise in subsonic compressible flow, where the infection 
spreads in some cases so that the formal second-order solution 
is incorrect not only near the edges but over the entire airfoil 
surface. Thus, using the particular integral of reference 14, 
Harder and Klunker gave an expression for the second-order 
solution for any symmetric airfoil at zero angle of attack 
(ref. 16). However, they noted that their expression does 
not apply to round-edged airfoils, for which it contains 
diveigent integrals. A more deceptive defect appears if 
their expression is applied to a sharp-edged airfoil such as a 
biconvex section; then the predicted surface speeds are 
finite (except near the edges) but incorrect everywhere by a 
term proportional to AP. This defect arises from the fact 
that near sharp edges the first-order somce distribution is 
not approximately the airfoil slope, as is assumed in thin- 
airfoil theory. The second-order solution involves the 
derivative of the source strength which, as indicated in 
figme 3, has infinite peaks that are missed by thin-airfoil 

Source 
strength 


00 CO 


Us 

derivative 


Thin -airfoil 
approximation 

Fiqubb 3. — Source strength for biconvex airfoil. 
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theory. It is enough to take account of this shortcoming in 
even the crudest fashion. Thus, if the region of integration 
is extended an infinitesimal distance beyond the edges to 
include the ptilses (Dirac delta functions) of the thin-airfoil 
approximation. Harder and Klunker’s expression yields a 
solution that is correct to second order except in the vicinity 
of the edges. 

Keune has discovered an alternative pai-ticular integral 
containing the stream function rather than the velocitj’’ 
potential, and so has obtained another expression for the 
second-order solution (ref. 15). Because the tangoncy 
condition is one degree smoother for the stream function than 
the velocity potential, his expression yields the correct 
result (except near stagnation edges) for shai’p-odged shapes. 
It fails, however, for round-edged shapes, so that liis solution 
for subsonic flow past an eUipse is incorrect everywhere. 

Both these expressions can bo manipulated by partial 
integration so as to be correct except near stagnation edges. 
However, the result is simply that obtained by applying 
the second-order compressibility rule to the expressions for 
second-order incompressible flow. Hence these more serious 
difficulties are of no further concern here. They do serve, 
however, to worn of the danger of false second-order solutions 
in more compheated problems. 

The role of edge corrections. — Thin-airfoil theory fails 
near stagnation edges because there the basic assumption of 
small disturbances is violated. It might be feared that uni- 
formly valid solutions coifld be foimd only by abandoning 
that assumption, which leads to such great mathematical 
simplification. Fortunately, however, Riegels and Light- 
hill have shown that for round edges in incompressible flow 
aU the resiilts of small-disturbance theory can be salvaged. 
They have given simple rules to be applied to the formal 
thin-airfoU solution that render it uniformly valid near the 
edge. 

For present purposes, corresponding rules are required for 
subsonic as weU as incompressible flows, for sharp os woU 
as round edges, and for cambered round edges, for which 
Lighthfll’s nfle is correct only to first order. It is beUeved 
that neither Riegels' nor Lighthill's technique can bo ex- 
tended to these cases. Instead, a different technique is used 
here, which is particularly suited to the study of edges. It 
consists in comparing the exact and thin-airfoil solutions for 
simple shapes that approximate the airfoil in the vicinity of 
the edge. This technique was first apphed in reference 20 
to the surface velocity on airfoil sections, tlnee-dimensional 
wings, and bodies of revolution. It is reproduced here 
insofar as it apphes to second-order theory for airfoil sec- 
tions, and is extended to treat sinface pressures as well os 
velocities. 

Bound edges in incompressible flow. — Most subsonic air- 
foil sections have finite leading-edge radius, and many are 
actually analytic (except at the training edge) . This moans 
that all derivatives are continuous, so that, Avith s the 
abscissa measiued from the edge into the airfoil, the upper 
and lower surfaces are described by 


y=±To8^^+C,s±T^^'^+C^+ . . . 


(14) 
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Here the T, and are coefficients for the tliickness and 
camber, respectively (fig. 1). Thus To deteimines the nose 
radius, Oi the initial camber, and so on. 

Thin-ah’foil theory breaks down close to the edge where s 
is of the order of the leading-edge radius p. Now p is pro- 
portional to the square of the thickness ratio t, so that thin- 
airfoil theory faUs where s=0(t^). In this small region the 
airfoil is described to second order by the first two terms of 
equation (14). In terms of the leading-edge radius p and 
initial slope X of the camber fine, the airfoil is given by 

4 — i/2pj?-i-X.9 (15a) 

This is the equation of an incHned parabola, described in 
rotated coordinates by 

y= ±V2ps (15b) 


where the two origins of coordinates are separated by a 
negligible distance, and the difference between p and p can 
also be ignored. 



Figure 4. — Flow past inclined parabola. 


The exact velocity on this parabola in a uniform stream 
is found, from conformal mapping or otherwise, to be given 



where the signs refer to the upper and lower surfaces. Here 
a has the physical interpretation that the stagnation point 
lies at 3 = 0 ^ and ^=-a-^ (%• 4). When the flow past 
the parabola is related to that near an airfoil nose, a is some 
moderate multiple of the angle of attack measured from the 
"ideal” angle — the angle at which the stagnation point coin- 
cides with the vertex. The factor of proportionality depends 
on the entfie airfoil shape (and the trading-edge condition), 
but its value is not required here. Expanding this expression 
for small p/3 yields, to second order, the formal series 



a p 

vf ^ 


and, as it must be, this is the second-order thin-airfoil 
solution for a parabola. It is clear from this formal expan- 
sion how the singidaiity ai'ises at the leading edge. 

The ratio of the exact speed on the parabola to its formal 
series expansion serves as a multiplicative correction factor 
for othei- shapes having the same nose radius. Thus the 
second-order thin-airfod solution ga for any airfod of leading- 
edge radius p is. converted into an appro.ximation g that is 
uniformdy valid near the edge by 


45 


(18) 


Simplifying this insofar as possible Avithout destroying its 
validity near the edge, and retaining only second-order terms, 
gives the ride 


?_ 



(19) 


It might be supposed that since the airfod nose was fitted 
to second order, this rule yields a solution that is imiformly 
valid to second order. Comparison with various exact 
solutions indicates that this is true, in the sense that the 
velocity disturbance and hence the pressme coefficient are 
correct to second order everywhere (except near the trading 
edge, where additional modification is required). This will 
be indicated by replacing g by gs. 

The oblique abscissa 5 can be expressed in terms of the 
original abscissa s, since on the surface of the parabola 

5=siXV2p8-l- . . . (20) 


Hence the rule becomes finally 


gz _/ 8±XV^ 

U \8±\^/^+p/2j 


( 21 ) 


In the special case X=0 this reduces to Lighthdl’s rule 
(ref. 5). 

The corresponding rule for siuface pressure coefficient is 
found by proceeding in the same way with the exact pressure 
on the parabola, which is found from equation (16) using 
Bernoulli’s equation. Again expanding formally for small 
p/s and taking the ratio as a multiplicative correction factor 
yields the rule 


s±X-\^^-fp/2 


( 22 ) 


Airfods of the NACA four- and five-digit series are not 
analytic at the nose. Their thickness distribution T{x) 
consists initially of the ordinates of a parabola minus those 
of a wedge, so that the airfoil is described by 


y= ± ToS^+CisT Txls|-f. . . (23) 


(17) 


These airfoils are fitted only to first order by an inclined para- 
bola, and it follows that the preceding rules render the thin- 
airfod solution uniformly valid only to first order near the 
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edge, though it remains correct to second order elsewhere. 
An appropriate second-order rule would require finding the 
conformal mapping for a parabola minus a wedge. 

First-order theory; Eiegels’ rules. — In. the first-order 
theory, the terms in X can be neglected, so that the rule for 
pressure, equation (22), simplifies to 


/O — ° ^ 


(24) 


The last term in equation (21) is also of second order, so that 
the rule for speed becomes 


U~(s+p/2) U 

This is not precisely Riegels’ nile. However, for a parabola 

where ri is the angle of the surface. Hence an alternative 
form of equation (25) is 




(27) 


and this is Riegels’ rule (ref. 2). The corresponding rule for 
pressure is 

Op=(cosu)»(7^ (28) 


These alternative forms suffer the defect of falsely implying 
that the nose exerts an influence even at remote points if 
the local airfoil slope is appreciable (for example, at the 
trailing edge). To this extent they fail to render the solution 
uniformly valid. However, it happens that they are exact 
for elhpses as well as parabolas, and are accordingly much 
more accurate for most airfoils, as indicated by the example 
of figure 5. 



Figubb 5. — CompariBon of Lighthill's and Riegels’ rules for incom- 
pressible flow past NACA 0012 airfoil. 


The second-order rules can be manipulated into the form 
of Riegels’ rule (ref.20), but the sUght advantages of accuracy 
and simphcity then scarcely offset their defects. The same 
is true of the rules for compressible flow given later. Hence 
Riegels’ rules are recommended only for first-order incom- 
pressible-flow theory. 

Round edges in subsonic flow. — The previous rules can 
be extended to subsonic speeds simply by considering sub- 
sonic rather than incompressible flow past an inclined 
parabola. Although no exact solution of this problem is 
known, existing approximate solutions by the Janzon- 
Rayleigh method are probably sufficiently accurate, at least 
at the ideal angle of attack, and could, in principle be refined 
indefinitely. Alternatively, one could use experimental 
measurements on a parabola. 

The problem is defined by p and a and the free-stream 
Mach number M. (The adiabatic exponent 7 also enters, 
but is assumed fixed at 7/5.) Therefore, dim ensional reason- 
ing shows that the surface speed is given by 



where, as usual, the upper and lower signs apply to the upper 
and lower surfaces. (Choosing p/2 rather than p as the ref- 
erence length leads to later simplification.) Expanding this 
function formally for x large compared ■with p and a’ would 
yield the thin-airfoil series, which is, to second order 

( 30 ) 


where Ki and K% are the compressibihty factors of equation 
(12c). Again the ratio of Q, to its series expansion servos as 
a multiphcative correction factor to be applied to any round- 
nosed airfoil. Simplifying as before, and replacing if by s 
according to equation (20), gives the rule* 


q% n/«±xV^ , a 





i 8 



Ki-1 

2 



(31) 


where qi/U is the first-order thin-airfoil solution. Hero a 
must be identified as the coefficient of ±Kif-y}s in the second- 
order solution, as is clear from equation (30). (The physical 
interpretation of a sho'wn in fig. 4 is valid only at zero Mach 
munber.) 

The corresponding rule for pressiue is foimd by consider- 
ing the exact pressure coefficient for the parabola, wliich 
must have the form 


> As Af— *0 this reduces not to the rale given prevlou 8 l 3 r /or Incompressible /lotr, but to on 
altematiye that b equivalent up to terms of second order. The difference arises trom tho 
laot that in the Incompressible case the dependence of g/U upon a Is given explicitly by tho 
factor (ld=a/V^) In equation (16), which Is used to cancel a corresponding term In tho 
denominator of equation (18). 
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Its series expansion would give the thin-airfoil result 

(33) 

Hence the rule for any airfoil is 


where 



(34a) 


(34b) 


In this case a is the coefficient of ^2£^i/Vs in the second- 
order pressm-e coefficient. 

Imai has recently calculated the Janzen-Rayleigh solution 
for a parabola at the ideal angle of attack including terms in 
Af* (ref. 31). Thus, his results give 


Q(^,Q,M'^=g_o+M%+M%+0{M<^) (35) 


whore the g„ are increasingly complicated fvmctions of 
ff/(p/2). In reference 20 an attempt was made to increase 
the accuracy of this approximation by modifying it so that 
for large 3f/(p/2) it tends exactly to the second-order thin- 
airfoil solution of equation (30). However, the third-order 
tliin-airfoil solution has since been calculated by Kaplan 
(ref. 32): 

g-l X, (x)+I K, [(0 + 

where 



r T-kiiip 

f T+l ilPN 

\ (8-AP) 

2 /3® 

‘ 8 /3® ' 



(It is a matter of taste whether the logarithmic term in- 
cluded here is regarded as being of third or fom1;h order.) 
The comparison with this result shown in figure 6 suggests 
that the modification was detrimental, and that Imai’s 
solution is adequate for practical purposes.® A short table 
of the function Q[3/(p/2), 0, M\ calculated from equation (35) 
is given below: 


M 

B 

B 

B 

2 

B 

10 

20 

60 

0 

0 

a 6774 

0.7071 

a 8165 

aS044 

a9636 

a 9769 

a 9001 

.4 

0 

.6670 

.6874 

.8002 

.8829 

.9473 

.9724 

.9886 

.6 

0 

.5449 

.6762 

.7808 

.8763 

.9431 

.9690 

.0874 

.6 

0 

.5205 

.6502 

.7768 

.8848 

.9371 

.0664 

.0358 

.7 

0 

.6104 

.6388 

.7674 

.8508 

.9200 

.9616 

.0835 

.8 

0 

.4870 

.6133 

.7338 

.8324 

.9182 

.0551 

.0805 

»S5 

0 

.4737 

.6982 

.7197 

.8213 

.9116 

.9512 

.0788 

.9 

0 

.4500 

.6816 

.7030 

.8088 

.6042 

.9466 

.9765 



The corresponding table of the fimction P[?/(p/2), 0, 
appearing in equation (34) is: 


P 



M 

0.5 

1 

2 

4 

10 

20 

50 

0 

0.3333 

a6000 

a6667 

asooo 

0.9001 

a 0524 

aoso4 

.4 

.2024 

.4444 

.6018 

.7339 

.8502 

.9016 

.9402 

.6 

.2823 

.4010 

.6500 

.6779 

.7958 

.8507 

.8050 

.6 

.2180 

.3391 

.4703 

.5881 

.7026 

.7687 

.8073 

,7 

.1607 

.2523 

.3560 

.4631 

.6523 

.6044 

.6517 

.8 

.0018 

.1466 

.2111 

.2742 

.3427 

.3797 

.4165 

.85 

.0578 

.0932 

.1367 

.1783 

.2256 

.2518 

.2774 

.9 

.0284 

.0463 

.0682 

.0906 

.1162 

.1305 

.1437 


For other angles of attack, the function Q to order 
can be extracted by a limiting process from Kaplan’s Janzen- 
Rayleigh solution (ref. 33) for an inclined ellipse. This 
gives, with ?/(p/2_)=A', al-\Jp]2=A 


1+^T/ (V-- 
1+A 


[(VA:-biAX-iA)ln 


1+x 


(l — tan 



1 Tho tabic on p. 11 of ret 20 also Involves an error In computing the of eq. (36). 


(36) 
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The corresponding approximation for the pressure coefficient, is 


J1{X, A, ]>,£)■ 


l-2A-ylX-A^ , 
1+X ^ 




i (l -2A VX-A*) =+(l -A^) {X+A-yfT} -A{X+A^) (V^+A) + 


(1+A*)(VX+^) 

l+X 


[(V^+|aZ-|^) In i^+(i-X+2^V^ tan-‘V^] 


(37a) 


and for the fimction P, according to equation (34b) 

n^XyAM) 


P(X,A^i)= 




(37b) 


The error involved in retaining only terms of order 
can be estimated from figure 6 in the special case a=0. At 
other angles of attack the error may be greater; in particular, 
neither the Janzen-Rayleigh expansion nor the thin-airfoil 
expansion is believed to converge if the local Mach number 
exceeds unity. 

Sharp edges. — The coiTections for sharp edges can be 
found by considering flow in an angle. At a trailing edge 
mth Kutta condition enforced, or a leading edge at ideal 
angle of attack, the surface and the dividing streandine meet 
at slightly less than a straight angle, as indicated on the 




'~^s^-,^rnr77777777777777 

I 



Figubb 7. — Flow near sharp edges. 


right of figure 7. For incompressible flow, the surface speed 
is found from conformal mapping to be given by 

(38) 

where 5 is the semivertex angle, and s is again measured into 
the edge. In fixin g the constant c, the difficulty that in the 
angle flow the velocity increases indefinitely upstream can be 
circumvented by requiring that at any point the velocity 
must approach that of the free stream as the angle 5 tends to 
zero. Thus it is seen that c is imity except for terms of 
order 5. 

The connection with thin-airfoU theory follows from the 
fact that for small e 

s‘=l+e In s+i log*s+. . . (39) 


though this relation is not muformly valid near s=0. 
Hence, the thin-airfoil series for flow in an angle is, to second 
order. 

(40) 

and it is clear how the spurious logarithmic singularities arise. 

Comparing these expressions gives a rule that renders the 
second-order thin-airfoU solution for any sharp-nosed profile 
uniformly valid: 

a_,A[g-ifiln.-5(ln,-lln',)] (41) 

For a leading edge not at the ideal angle of attack, the 
flow includes a circulatory component, as indicated at the 
left of figure 7. For incompressible flow, conformal mapping 
gives the surface speed associated with this component os 

^=±ttS (42) 

and its thin-airfoU expansion is 

Comparing these yields a rule for correcting the ch’culatory 
component: 

a=,S^(§-|;gta.) (44) 

The second-order thin-airfoU solution can be treated by 
splitting off the terms that are singular at least as near 

the leading edge, applying equation (44) to this circulatoiy 
component, applying equation (41) to the remainder, and 
recombining. 

These rules could be extended to subsonic speeds, in the 
case of ideal angle of attack, by calciUating the Janzen- 
Rayleigh solution for flow in an angle. However, at other 
angles the Janzen-Rayleigh approximation certainly fails, 
because it would predict infinite speeds that are tolerable 
only in an incompressible fluid. In any case, the correction 
is negligible for most practical purposes, because it is signifi- 
cant in only a minute neighborhood of the edge, and, fur- 
thermore, sharp edges are usually trailing edges, in wliich 
case the detaUs of the flow are masked by viscous effects. 
For these reasons, no correction for sharp edges is included 
in the computing scheme given later. 

Combined edges. — AirfoUs with two stagnation edges are 
treated by applying the appropriate corrections in tmn at 
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each, edge by shifting the origin. Thus, consider an airfoil 
with round leading and trailing edges of radii pa and p* 
located at x=Xa and x=Xb, respectively, and initial and final 
camber angles and Assume that the Kutta condition 


is imposed at the trailing edge. Applying equation (31) 
twice, identifying s successively with x—Xa and Xt,—x, and 
then simplifying to keep no more than second-order terms 
gives 


Si— 

U~ 


X Xaziz\i-\/2pf,(x x^ , a Tir n ^6 ^ ^ nr 

® L v^' 


Qi -T- T7 ® 2l I -^ 2 / Pa I Pb \ I f TT i -^2 1\ o* 


(45) 


(The simpler form of this equation for incompressible flow 
is given as equation (24) of reference 20.) The correspon ding 
rule for pressiue coefficient is 


valid second approximation by applying equation (21) 
twice in succession, or using the combined rule of equation 
(24) of reference 20, which gives 


<^'P2=^[- 


x—Xa±'Knl2pJ^x—x^ , a 

\ 


X 


Pb!^ J 


Similar rules can be foimd for combinations of a rormd 
and a sharp edge, or two sharp edges. For example, for 
incompressible flow past two sharp edges of equal angle, 
both Avith Kutta condition imposed (as for a symmetrical 
biconvex airfoil at zero angle of attack) and located at 2 = ± 1 , 
the combined rule has the form of equation (41) with « 
replaced by (1—®*). 


EXAMPLES: COMPARISON WITH EXPERIMENT AND OTHER 

THEORIES 


INCOMPEESSIBI/E FLOW 

It has been seen that the solution for subsonic flow depends 
on that for incompressible flow. It is therefore pertinent to 
test the second-order theory in the case of incompressible 
flow, where it can be checked against the exact results of 
conformal mapping. 

Ellipse. — Consider an ellipse of thickness ratio r with the 
interval — l<a;<l as chord line. It is described by 

J/=±r■^/^=^, -l<a:<l (47) 

Suppose that the Kutta condition is satisfled — the rear 
stagnation point coincides with the end of the major axis. 
Then the flrst-order solution for surface speed is found, 
from equations (3), (4), and (5), together with Appendix 
B, to be 

g=l+r±a^i^ (48a) 


Proceeding with equations (6) to (9) gives the formal second- 
order result 




This can be checked by expanding the exact result, which is 

U 

The formal second-order solution clearly breaks down near 
the ends of the ellipse. It is converted into a uniformly 





Figueb 8. — Speed on 18-percent-thick ellipse at zero angle of attack in 
incompressible flow. 

These approximations are compared in flgure 8 with the 
exact solution for an 18-percent-thick ellipse (which has 
nearly the same nose radius as an NACA 0012 airfoU) at 
zero angle of attack. The precipitate descent of the formal 
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second-order solution toward negative infinity is just dis- 
cernible near the nose and is eliminated by the edge cor- 
rection. 

Symmetrical JoukowsM airfoil. — ^To second as weU as first 
order a symmetrical Joukowski airfoil of thickness ratio t 
is described by 


Ti=-^ t=0.7698t 
3^/3 


(50) 


By the foregoing procedure, the formal second-order solution 
is found to be 


(1 1 -f ? r [2-(l-fx) In (1 -l-a;)- (1 -a:) In (1 -x)] -h 
- 0rJ[3-^(l-a^-3(H-x)ln(l-t-a;)-3(l-a:)ln(l-a;)-|- 

i (1 +3x) (1 -hx) In* (1 -hr) -i-i (1 -3a:) (1 -x) In* (1 -a;) -j- 
|(1 -r*)ln(l -f a:)ln (1 -a;)J_ia*]-±o^^/S(i +a:)^j^l - 
^ [2(l-fa;)ln(H-a;)-(l-t-2a;)ln(l-a:)-4] j- (63b) 




^(l-f2r)*T 


where the first three terms give the first-order solution. 
Correcting this by means of equation (21) with 8=l-l-x 
and p=4ti* (and X=0) gives the miifonnly vahd second 
approximation 

I- VnSk- V(l-2x)'± 

In figure 9 these approximations are compared with the 
exact solution (ref. 34) for a 12-percent-thick section at 
zero angle of attack. The effect of the edge correction on 
the second-order result is not discernible to this scale. 

Biconvex airfoil. — To second order a symmetrical biconvex 
airfoil of thickness ratio t bounded by either circular or 
parabolic arcs is described by 

y=±r(l-x*), -l<x<l (52) 

The formal second-order solution is found to be 




(53a) 

In deducing from this a uniformly vahd second approxi- 
mation, the terms independent of a are treated by the rule 
for combined equal sharp edges that was described just 
after equation (46), with 5=2 t. The terms in a are modi- 
fied according to equation (44) with s=l-|-x. The result is 



Figueb 9. — Speed on 12-percent-thick symmetrical Joukowski airfoil 
at zero angle of attack in incompressible flow. 
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FiauRE 10. — Speed on 18-percent thick biconvex airfoil at zero angle of 
attack in incompressible flow. 


These approximations are compared in figwe 10 with the 
exact solution (ref. 36) for a circular-arc airfoil 18 percent 
thick at zero angle of attack. Although the vertex angles 
are large in this example, the edge correction is appreciable 
in such a small neighborhood of the edge that it would be 
invisible even on a much larger plot. 

NACA OOXX airfoils. — Symmetrical airfoils of the NACA 
OOXX family (such as the NACA 0012) are naturally de- 
fined for the interval 0<a:<l. The airfoil of thickness 
ratio T is described by (ref. 36) 

0<a:<l (54) 

where 

6i= 1. 48450 

6j=-0. 63000 
64= — 1.75800 
^6= 1. 42160 

6s=— 0. 60750 
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With the aid of Appendix B the first-order solution is foimd 
to be 



(366-f268)a:-4b8K"J ±cc.yj^ (55a) 

in agreement with the result given by Goldstein (ref. 28). 
Applying Kiegels’ rule (eq. (27)) renders this a uniformly 
valid first approximation except" very near the trailing 
edge. 

The second-order terms in thickness, in addition to being 
very comphcated, involve integrals that apparently cannot 
be evaluated in terms of tabulated functions. Accordingly, 
the second-order terms have been calculated using the 
Kiegels-Germain numerical method discussed in Appendix 
C, with N=16. The accuracy of this approximation is 
assured by the fact that cruder approximations modify the 
numerical results only shghtly, as "will be seen in a later 
example. 

The formal second-order solution for surface speed there- 
fore has the form 

^=l+TQT±“.y/^^+T‘'QT-r±TQ:Q„— io* (55b) 

where values of Qr from equation (55a) and apprommate 
niunerical values of Q.„ and Q.ra are 


X 

Or 

Orr 

Or, 

X 

Or 

Orr 

Ora 

a025 

1.043 

“9.00 

8.80 

aso 

0.900 

-0.135 

as2 

.05 

1.838 

“3.35 

5.55 

.60 

.697 

-.220 

.11 

.10 

1.714 

“1.00 

3.25 

.70 

.485 

-.315 

-.08 

.20 

1.510 

“.090 

1.65 

.80 

.238 

-.410 

-.23 

.30 

1.309 

.010 

LOO 

.90 

-.m 

-.420 

-.33 

.40 

L106 

-.060 

.58 

.95 

-.440 

-.380 

-.34 


Applying equation (21) "with p=l. 10187 t*, 5=1.16925 t 
(and X=0) yields a unif ormly valid approximation. How- 
ever, as disciissed previously, the curvature of the profile 
does not vary continuoiisly near its nose, so the result is 
only a first approximation there, though a second approxi- 
mation elsewhere. 

The various approximations are compared in figure 11 
■with the resiilt of a “long and elaborate calculation” by 
conformal mapping for the NACA 0012 airfoil that is given 
by Goldstein (ref. 37). Again the effect of modifying the 
second-order solution is indiscernible. Also sho"wn is the 
“exact” solution tabulated in reference 38. The agreement 
between the first-order solution "with Riegels’ rule, the 
second-order solution, and Goldstein's calculation leaves 
little doubt that his is the more accurate of the two “exact” 
solutions. 

COMPRESSIBLE FLOW 

When extended to subsonic compressible flow, the preced- 
ing examples can all be compared "with other theories or -with 
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theory, Riegels' rule 
t-order theory 
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— 

LJ 


Fraction of chord 


Figubb II. — Speed on NACA 0012 airfoU at zero angle of attack in 
incompressible flow. 


experiment. As before, the comparisons ■will, for sunplicity, 
be made only for zero angle of attack. 

Ellipse. — ^Applying the second-order compressibility rule 
of equations (13) to the incompressible solution of equation 
(48b) gives as the formal second-order solution for the speed 
on an elliptic cylinder 

g-i+jr. (r±»7^+i '■ ^“0=^ 

For zero angle of attack the maximum speed, occurring at 
midchord, is given by 

§)_^=l+Arir-f|(.B:.-l)r» (57) 


in agreement -with the result of Hantzsche and Wendt (ref. 7) 
The corresponding third-order result has been given by 
Hantzsche (ref. 8); his expression can be simplified to 


»■)}'’ -i['+ 

f (l+l)(8-^')] ;+0(^^) (58a) 


where 


T-l-1 

2 


(68b) 


The terms in r^ln r are found correct together with those in 
T® in the second iteration; it is a matter of taste whether they 
are regarded as being of third or fourth order. The solution 
to T® has subsequently been given by Hasimoto (ref. 39), and 
although typographical errors unfortimately appear in his 
equations, his nmnerical results agree wth those calculated 
from equation (58a). 

Values of the maximum speed ratio calculated from these 
and other approximations for a 10-percent-thick ellipse at 
zero angle of attack (whose critical Mach munbor is about 
0.80) are: 



AI-O.TO 

iU=0.76 

A/- 0.80 

Flret-order theory (or Prandtl-QlauertnUe applied 

1.140 

1.151 

1.107 


1.149 

1.163 

1.184 


1.148 

1.163 

1.185 

Thlid^jrder theoryr 

1.150 

1.167 

1.190 


1.149 

1.166 

1.194 



(Here the Kdrmdn-Tsien rule has been applied to the pressure 
coefficient for incompressible flow, and the speed ratio then 
calculated from Bernoulli’s equation.) It might have been 
anticipated that the second-order theory is more accurate 
than any of the compressibflity correction formiflas such as 
the Kdrmdn-Tsien rule, because it allo\vs for a dependence 
on the particular airfoil shape and on the value of 7. How- 
ever, in this example the results of second-order theory 
and the KArmfin-Tsien rule are practically identical. 

In the same way the second-order solutions are readily 
calculated for the Joukowski and the biconvex airfoils, and 
are foimd to agree -with the results that Hantzsche and 
Wendt obtained by laborious analysis. 

NACA 0012 airfoil. — The formal first- and second-order 
solutions for NACA 00 XX airfoils in subsonic flow are 
easily obtained from equations (13) and (55). The second- 
order solution can then be rendered uniformly valid near the 
nose using equation (31), although again the modification is 
significant in only a very small region of the nose. 

For the NACA 0012 airfoil at zero angle of attack, Em- 
mons has calculated the flow field at Mach munbers of 0, 
0.70, and 0.75 using the numerical relaxation method (ref. 
40). The last of these Mach numbers is supercritical, so 
that the flow contains shock waves, and is beyond the scope 
of the present theory. The pressure distribution calculated 
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Fiqueb 12. — Pressure coefficient on NACA 0012 airfoil at zero angle of 
attack in subsonic flow. 


by the relaxation method for M=0.70 is compared in figure 
12 with the results of first- and second-order theory and 
various other approximations. The relaxation solution for 
incompressible fiow is also shown in comparison with Gold- 
stein’s “exact” solution, and is seen to be inaccurate near the 
nose. The solution for M=0.70 probably contains similar 
inaccuracies; however, just as for the ellipse the pressure 
coefficients calculated by second-order theory may be shghtly 
less negative than the true values near the minimum. 

Experiments on NACA 0016 airfoil. — ^Experimental pres- 
sure distributions in two-dimensional fiow over the NACA 
0016 airfoil at high subsonic speeds are reported in reference 
41. For zero angle of attack, the critical Mach mmaber is 
approximately 0.70. The measurements at this Mach 
munber are compared in figure 13 with the results of first- 
and second-order theory and of the two common compressi- 
bility correction formulas applied to the incompressible 
flow values tabulated in reference 38. Unfortunately, 
the model was imperfectly constructed, and the ordinates 
were inaccurate near the nose and midchord. Otherwise, 
the measured pressures are in satisfactory accord with 
either second-order theory or the results of the KArmkn-Tsien 
rule. 

Tomotika and Tamada’s airfoil. — Using the hodograph 
method, Tomotika and Tamada have calculated the flow 
past a certain family of symmetric airfoils (ref. 42). As 
usual in hodograph solutions, the airfoil shape varies some- 
what with free-stream Mach number. 'The critical Mach 
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Figurb 13. — Comparison of theoretioal and experimental pressure 
distributions on NACA 0015 airfoil at M=0.70, zero angle of 
attack. 


number is 0.717, and the corresponding shape is shown in 
figure 14 together with the surface speeds predicted by 
various theories. 

For mathematical simplicity, Tomotika and Tamada have 
adopted a hypothetical gas, which is fitted at Mach numbers 
zero and unity to a polytropic gas having 7 = 7 / 5 . At any 
intermediate Mach number, however, the hypothetical gas 
corresponds to a polytropic gas whose 7 is greater than 7/5, 
reaching a maximum value of 1.91 at M=0.78. To second 
order, any such hypothetical gas is equivalent to a polytropic 
gas having the vedue of 7 corresponding to the free-stream 
flow, given by 



where p is the density and c the speed of sound in the hypo- 
thetical gas. For Tomotika and Tamada’s gas with M= 
0.717, that value is 1.82. Actually, the second-order solution 
depends so slightly upon the value of 7 that the change from 
7 = 7/5 to 7 = 1.82 increases the maximum value of gJUhj only 
two parts in a thousand. However, the nonpolytropic 
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FiaiTHB 14. — Speed on Tomotlka-Tamada airfoil at Af=0.717. 


foil at any angle of attack. It requires a knowledge only of 
the airfoil ordinates at 7 (or 15) points along the chord. It is 
based on the foregoing theory together with the numerical 
method of Riegels and Germain that is discussed in Appendix 
C. It includes the correction for a round leading edge, but 
not for sharp edges. 


COMPUTING PROCEDURE WITH iV=8 (OR 16) 

(1) Tabulate the ordinates and Fi of the upper and 
lower surfaces at the 7 (15) pivotal points x„ hsted in table I. 
(The X axis must pass through the leading and trailing edges 
with the leading edge at x=0 and trailing edge at a:=l.) 

(2) Calculate the corresponding values of 

T=i(r„-F,), C=i(T„+r,) (59) 


(3) Using the influence coefficients of tables II, III, and IV, 
calculate 


7(16) 

ni * 7(15) 

U 

a-1 


C'=^UC. - 

j-i 


7(16) 

S 9n,T, 
*•1 

706) 

i-1 


(60) 


nature of the hypothetical gas must be considered in con- 
verting Tomotika and Tamada’s values of stuface speed 
(which are referred to the critical speed) to the form g/U. 
At the free-stream Mach number of 0.717, the ratio of 
freestream to critical speeds is 0.769 for the hypothetical 
gas, compared with 0.748 for a polytropic gas having 7=7/5 
(and 0.774 for a polyfropic gas having 7=1.82).* Thus, as 
Tomotika and Tamada emphasize, their results should be 
regarded as exact for the hypothetical gas; and the present 
method should he regarded as giving the true second approxi- 
mation for the same gas (which to second order is equivalent 
to a polytropic gas with 7=1.82). 

The agreement between the hodograph method and 
se<xmd-order theory is satisfactory. Third- and higher-order 
terms are seen to increase the surface speed near its maximum 
and reduce it slightly elsewhere, which is the efltect of 
third-order terms calculated by Asaka for a biconvex airfoil 
(ref. 43). 

As shown in figvue 14, the residts of second-order theory 
and the Kfirm4n-Tsien rule nearly coincide except over the 
forward portion, where the second-order theory is apparently 
more accurate. The second-order solution calculated with 
the present computing scheme is in close agreement with an 
impublished second-order solution carried out by Naruse. 

PRACTICAL NUMERICAL COMPUTATION 


(4) Using table I (with a in radians), calculate 


(5) Calculate 


the '^1* I jl—x 

U~ U 






rp rji t '^le pt 

^ -2—^ - I Jj 


(61) 


(62) 


U a=l ^ U 


(6) Using the influence coefficients of table II, calculate 

(63) 

(7) Using the compressibihty factors of table V, calculate 

(^Jf.CC'+TT"+\ O't+l T'^^+ - 

^ Kv)‘+(v)>[^' TT+^" (?'+«’"+ 


or 


C''T+C'T'^+{K2-1) ^ 
C^=-2K, [2 ^±2 


(64a) 


The following computing procedure yields the second-order 
subsonic solution for the surface speed or pressure on any air- 


2{C±T){0''±T'')+{C'±T'f^(^±^ (64b) 


< The author Is Indebted to Mr. H. Naruse of Tokyo Metropolitan TTniversIty for having 
pointed out these Ideas, some of which stem In turn from his coUeegues, L Imal and S. 
TaknmL 


The ± signs refer to the upper and lower surfaces of the 
airfoil. 
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Correction for round leading edge. — the airfoil has a 
round leading edge, continue as follows: 

(8) Using the influence coefficients of table I, calculate 


V 2 


7(161 

ai*=^d„Cn, 

n»l 

(9) Calculate 


n=l 

7(16) 

G^*= ^ f AnC^n 


(65) 

( 66 ) 


®=V2(“+ai*+(i2*)+2 (67) 


( 68 ) 

(10) At each pivotal point on the upper and lower surfaces 
calculate 

(69) 

and then calculate 

Q(Xn,dcA,M) or P(X^,:kA,M) (70) 

from equation (36) or (37) 

(11) Using the compressibility factors of table V, calculate 
at each point 



U' 



Ki-1 

2 



(71a) 

or 

C\=PC^^ (71h) 


EEMAEKB 

1. The summations of steps (3), (6;, and (8) are conveni- 
ently carried out by tabulating C and T in columns that can 
bo matched with successive columns of tables II and III 
while cumulative mtfltiplication is carried out on a desk 
calculating machine. 

2. If the airfoil slope and second derivative can be found 
more directly, step (3) can be simplified by omitting the 
calculation of T', O', T" , and 0" ■ 

3. If the leading-edge radius p and initial camber angle X 
are known, omit their calculation from step (8) . 


4. Seven pivotal points yield sufficient accuracy for most 
purposes. If conditions near the nose of a thin airfoil are 
of interest it may be necessary to repeat the computation 
using 16 pivotal points. 

5. The above scheme is designed for calculating a single 
case. If the same airfoil is to be calciflated at more than two 
angles of attack, it is economical to subdivide the computing 
scheme to separate terms in a and «*. Similarly, the scheme 
should be subdivided if more than two thickness or camber 
ratios are to be calculated for the same family of airfoils. 

6. Por NACA airfoils T is the basic thickness and 0 the 
camber line. To second order it is immaterial that the thick- 
ness is added normal to the camber line rather than to the 
chord line. 

7. In step (10), if A=0 the values of Q orP can he taken 
from the tables following equation (35). 

EXAMPLE 

The following table gives the complete computing form 
for calculating the first- and second-order increments in 
surface speed for an NACA 00 X X airfoil (of imit thickness 
ratio) at zero angle of attack and zero Mach number: 


n 

X 

T 

«U A^i 

U U 

•p 

T" 

T, 

Uu 

U 

ATj 

V 

1 

a03806 

a 26316 

L8938 

3.0914 

-54.721 

a40S32 

4.8(M0 

-4.8180 

2 

.M6I6 

.44238 

L6168 

.8662 

-a 628 

.71478 

a 1691 

-.4062 

3 

.30366 

.49690 

L2906 

-.0236 

-a 636 

.W617 

L8200 

.0028 

4 

.50000 

.44051 

.9053 

-.«46 

-L989 

.39379 

.5H9 

—.1840 

6 

.09134 

.30843 

.4909 

-.8063 

-1,(H8 

.16142 

-.3081 

-.3063 

6 

.86366 

.16199 

.0824 

-L0163 

-L674 

.01336 

-. 6913 

-.4293 

7 

.901^ 

.04499 

-.6951 

-L0S79 

L643 

-.02677 

-.9388 

-.3230 


The accuracy of this solution with only 7 pivotal points is 
indicated by comparison with the following values, which 
were obtained analytically for Agj/U and with 15 pivotal 
points for AgJU: 


X 

ail 



U V 

U 

a 14645 

1.6166 

-a 4009 

.50000 

.0003 

-. 1348 

.85355 

.0726 

-.4239 


It is seen that the solution using 7 pivotal points yields ample 
accuracy. 

Ames Aebonauticaii Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., March 12, 1956 
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A 


bn 

C(x) 

CM 

Op 

^Cp, 

d'n»j ^nxj 
fns} Qntj ^n* 

dn 

M 

fn 

I 

Ki 

IC2 


h 


M 

N 

P 


Q 

Qt> Qtt> Qra 
2 

A2i, A2a 

R 


T(x) 

TM 


APPENDIX A 


PRINCIPAL SYMBOLS 


al4^ 

factor proportional to angle of attack meas- 
ured from ideal angle 

coefficient of x”/* in series for airfoU ord- 
inate 

camber of airfoil 

camber of fictitious airfoil in second-order 
solution 

surface-pressure coefficient 
first-order smface-pressnre coefficient 
second-order increment in surface-pressure 
coefficient 

iufluence coefficients for calculating velocity, 
slope, and second derivative of airfoil 
ordinate 

influence coefficient for calculating a 
analytic function of complex variable 
/o„, influence coefficient for calculating ini- 
tial camber angle 

imaginary part of /(a) on unit circle 


first-order compressibility factor, 




second-order compressibility factor, 
(T-M)Af*d-4j8" 

413^ 

coefficient in trigonometric polynomial ap- 
proximation to 0 
free-stream Mach number 
number of subdivisions of chord line in 
numerical integration 

ratio of pressiue coefficient on parabola in 
subsonic flow to second-order thin-airfoil 
value 

surface speed ratio on parabola in subsonic 
flow 

(Seeeq. (55b).) 

flow speed on surface of airfoil 
first- and second-order increments in 2 


real part of/(z) on imit circle 
influence coefficient for calculating edge 
radius 

abscissae measured from edge positive into 
airfoil 

distance from round edge measured along 
initial tangent to camber line 
thickness of airfoil 

thickness of fictitious airfoil in second-order 
solution 


U 

U 

u 

V 
X 

X 

Y(x) 

Tp,Y 

y 

z 

a 

/3 

^P 

y 

yp 

5 

e 

V 
X 

Mp 

n 

p 


Tl 

<p 

( )a 
( ). 
( )c 


( )if 
( )» 

( ). 

( )p 

( ). 

( )« 

( )o 

( )l 
( )* 

( )' 
(-) 


coefficient in trigonometric polynomial ap- 
proximation to T 
free-stream speed 

velocity pertm-bation parallel to chord lino 
velocity perturbation normal to chord lino 
abscissa 
s/(p/2) 

ordinate of airfoil 

ordmates of upper and lower surfaces of 
airfoil, respectively 
ordinate 

complex variable 
angle of attack 

coefficient in munerical calculation of u 

adiabatic exponent of gas 

coefficient in numerical calculation of Y' 

semivertex angle of sharp edge 

polar angle 

angle of airfoil surface to chord lino 
terminal angle of camber fine to chord lino 
coefficient in numerical calculation of Y" 
pressure coefficient on parabola in subsonic 
fiow 

radius of round edge 
airfoil thickness ratio 
4 

— pT for Joukowski airfoil 

3V3 

pertm-bation velocity potential 
perturbation stream fimction 
value at leading edge 
value at trailing edge 

component associated ivith camber and 
angle of attack 
value at Mach number M 
value at mth pivotal point coimted from 
trailing edge 

value at Tith pivotal point counted from 
leading edge 

index of summation, counted from trailing 
edge 

index of summation, counted from loading 
edge 

component associated with thickness 
value at zero Mach munber 
first-order approximation 
second-order approximation 
derivative 

uniformly vahd approximation 
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APPENDIX B 

AIRFOIL INTEGRALS 

The follomng are the Cauchy principal values for x* < 1 : 


2 . 

5. £ 

J-i *— f J -1 x— ? n 

. £ ^ — ^=0 

• f 7=4^.^=-- 

• £vpfe)‘*--“ 

«=-<-+!) 

• "(*‘+l) 

“■ £vr=fczj, ■<«=-' K>^+l) 


10 


£ ri ( iin] 1(3) • • . (n 2) 

2 i 2(4:) .. . (71-1) 


14- £Cr*=-* 

47. £2S,,=,(,_4,_4) 
19. £^‘if-Ki»=i±j-'-) 


20 


1 1 1+^ 


4-1. £ *-| 14 (w^:-^)-2 


22 




2x 


23 


• £^^*‘'f-l*‘('“’Si-'*)-9l4- 






n 


on— 




Vi“H?(^ — f) Vi'h'^ -\/^''~Vi"l“® 

. (£' ^H±i Vi+^ ln ^-t->^ -2V2 

J_1 x-f ^ ^ Y2-VT+^ 


In 


1+f 


27. 


Vl— a:* 


. ^=i=i- cif= 

J-iVT=F(a;-f) 

Xa 

■ ^-iVi^Ca;-?)*^' 
tsinkhl 

29. f 

Xivr=F(a;-f) V Vi=^/ 


28 


Vi^ 


1+f 


30 


Vl— a^C*— f) V -^1—7?) 

31. 

J-iVl=i^(a:-f) V3 VI=^/ 


32 


. f 7^=^ r i?- 

vr=?(a;-f) J -1 VI=Ftx-f) 


vU-C-D-T 0-l) ! S - 


(-1)>(1)(3) ■ ■ ■ 
2'(2v+l)(v!)*0-l 


1— 2r 


(i) 
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APPENDIX C 


THE KIEGELS-GERMAIN METHOD 


The numerical procedure introduced by Riegels (ref. 2) 
and Germain (ref. 21) can be adapted to give approximately 
the thin-airfoil velocities on any profile in terms of its ordi- 
nates at certain fixed points. Thwaites has appbed Ger- 
main’s procedme to thin-airfoil theory (ref. 23), and Weber 
has systematized the calculation of the slope and surface 
velocity (refs. 24 and 25). Here we must also find the 
second derivative. In applying the edge correction at a 
round nose we also require the edge radius. It is convenient 
to derive all these results from Watson’s analysis (ref. 22). 

Let/( 2 ) be regular within the unit circle, and on the unit 
circle have the form 

J{e*^=R{e)+U{6) (Cl) 

(Our R and I are Watson’s ^ and e.) Then following Ger- 
main, Watson approximates to R by the trigonometric 
polynomial 

N-l 

R{6) COS rd-\-tr sin rd) cos Nd (C2) 

r-l 

which can be made to coincide with R at the 2N equally 
spaced pivotal points 6=B^=mwjN. Thus he derives 
approximate formulas for / (aside from a constant), R', F, 
SR, and SI hi terms of the values of R at the pivotal 
points times fixed influence coefficients. 

In thin-airfoil theory the complex perturbation potential 
tp+irp is regular outside the unit circle in the absence of circu- 
lation. Inversion shows that this involves a change in sign 
of either the real or imaginary part, since /(«“**) 

Hence {<p, — Tp) or (^,p) may be identified with {R,I). 

In thin-airfoil theory the tangency condition on the per- 
turbation stream function 4' is 4'= — where Y is the air- 
foil ordinate. Therefore, in order to obtain a solution in 
terms of the airfoil ordinates we identify (F,— p) with {R,I). 
It is assumed throughout that the x axis passes through the 
leading and trailing edges. 


STREAMWISE VELOCITY INCREMENT FOE a=0 

Let X run from 0 at the leading edge of the airfoil to 1 at 
the trailing edge, and 


x=- (1+cos 6) 


(C3) 


Then the streamwise perturbation velocity on the airfoil is 
given by 

u b<p b<p/Z)9 2 b<p .p . 

U~bx~dxlde~ sine be 

Now according to Watson’s equations (10), (24), and (27), 
in the absence of circiflation the values of bpfbO at the points 
e„ are given by 




1^’ 

0 , 


iV(l— cosSp) 


p ==0 

p=even, notO 
> p=odd 

( 05 ) 


Now since Pnf-p=Pp and Fo= F,=0, 

2N-1 N-l N-l 

^2Pp^m+p=^2Pp-vt7p'\'^2fip+mI'w-p ( 06 ) 

p-0 P-1 P=1 


Symmetrical airfoils. — ^For a symmetric airfoil F 2 «--»= 
— Y,= — T,. Then according to equations (04), (05), and 
(06) 

This form is convenient for calculating the c„p. It is also 
easily shown, using trigonometric identities, that 


r N 




sm e„ 

0 , 


p±m=0 

p±m=even, not 0 


( 08 ) 


‘at? a ^ p±m=odd 

N(cos e„—co3 BpY 


which is Weber’s result (ref. 24). Corresponding forms can, 
if desired, be found for the other influence coefficients dis- 
cussed later. 

Antisymmetric airfoils.— For a cambered airfoU of zero 
thickness Fw_p=Fp=C'p. Equation (06) gives 

2N-1 N-l 

Sflp2^«+p=S(flp-w+flp+n)Op (09) 

P-0 P-1 

This expression represents the velocity on the unit circle 
into which the airfoil is mapped. , The Kutta condition will 
be violated at the trailing edge of the airfoil unless the ex- 
pression happens to vanish for m=0. Adding a component 
of circulatory flow changes the velocity on the circle by a 
constant. Hence the Kutta condition is enforced by sub- 
tracting from the expression of equation (C9) its value at 
m=0, so that 


2N-1 N-l 

Zlflpr„+p^S(^P_„-l-/3p+„-2/3p)Cp (CIO) 

p-O p=l 

Hence, according to equations (C4) and (C5) 

(Lp=^^(P,.p,+Pp+„-2P,) (Cll) 
u/ p-o sm 

SLOPE OF AIRFOIL 

The airfoil slope is given by 


^ 

°° dx sin e de 


( 012 ) 


Now according to Watson’s equations (29), (31), and (34) 

0, fJ=0] 






y (013) 
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or, since 

dV\ -^-1 

yp-t<Xp — S ^p+nXiN-p (Cl4) 

Symmetric airfoils. — Using the symmetry conditions again 
gives for symmetric airfoils 

S empTp, 6mp= - (Tp-m+Tp+m) (Cl 5) 

Antisymmetric airfoils. — Similarly, for camber lines 

GU=^JmpOp, /»p=-i(7p+„-7,_«) (016) 

p=i sm (7* 


so that 

Setting m=N and 0 in equation (Cl4) and using the sym- 
metry properties of the thickness and camber distributions 
gives 

[e== -2 s y^-pT„ 2 S y^p (C24) 

V 2 P=1 P”1 

at the leading and trailing edges, respectively. 

SLOPE OF CAMBER LmE AT END 


SECOND DERIVATIVE 

The second derivative of the airfoil ordinate is given by 

4 /(P7 

/» ( 7/I» cot 6 Jy. 




dx^ sin* d\dd^ 


(017) 


An approximation for d?YldB^, which is required here, is 
found by extending Watson’s analysis for the first derivative, 
as he suggests. Following closely his section 2.4 gives, after 
some computation. 


W A"" S ^p^«+p> -j 


2 W *+1 

6 

-(-!)” 


1 — cos d 


, p=0 
■> Pt^O 


(018) 


Symmetric airfcils.' — ^By the foregoing procedure, it is 
found that for symmetric airfoils 


— S g„pYp 

p=i 


Smp- 


sin* 0, 


■ [(Mp-p,-Mp+J-COt 0„(7p-„+7p+J] (019) 


Antisymmetric airfoils. — Similary for camber lines 

N-l 

^pYp, 

p-i 


^p g [(A*p-in"l”Mp+jii) cotd„(yp^„ 7p+«)] (020) 


EDGE RADU 

If the airfoil has round leading and trailing edges, their 
radii are given by 


p= lim 2iW*, p= lim 2(l-ar)W* 

i-<0 i-*l 

Near the edges, to a first approximation 



( 021 ) 


( 022 ) 


1 Tbo expression given for in NAOA TN 3390 is incorrect, differing ffom that of the 
present eqimtlon (018) by l/2(— l)*2V". The error cancels out In the Ompf but the values of 
hmp tabulated there are affected* However, the correction has little effect (except to simplify 
the tables) slnco the approxlmatlou merely becomes exact for the first $ Fourier components 
rather than the first 7 (In the case 


Equation (016) is indeterminate at the ends of the airfoil. 
Applying L’Hospital’s rule to equation (012) gives at 8=0 

Yo'=-2^ (025) 

Then using equation (018) and the fact that YiN-p=Yp=Op 
for a camber line gives 

-S fpGp, fp ^fop= -4(-l)^ -- ^ ^ (026) 

1 1 — cos Op 

ANGLE-OF-ATTACK PARAMETER a 

The parameter a is the coefficient of the square-root 
leading-edge singularity iu the surface speed. There is a 
first-order contribution from u^JU (eq. (4)), and second-order 
contributions from UtJU (eq. (8)) and (0±T) (0" ±T”)-\- 
1/2(0' ±T')* (eq. (9)). 

For zero angle of attack, equations (04), (05), and (OlO) 
give 

Near the leading edge 

sin 0~ (028) 

so that (since Pp-n=Pp+n=^n-p) 

g (P^-p-Pp)Cp (029) 

Including the contribution of angle of attack from equation 
(4) gives for the first-order approximation to a 

tti=V2 dp=2(fi^_p — ^p) (030) 

For a round nose described by equation (15a), the terms 
(0±T) (0"±!r")-l-l/2(C"±T')* of equation (9) contribute. 



Hence the second-order increment in a is 

Aoa=V2 dLpG^p~\i) k- /fi (031) 

p=i V 2 
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TABLES 

The various influence coefl&cients have been calculated for 
iV=8 and A'^=16. The values have been checked by apply- 
ing them to a number of simple shapes — ellipse, parabolic 
arc camber line, etc. — ^for which the approximation of the 
airfoil by a trigonometric polynomial is exact. They were 
also checked against Weber’s values (refs. 24 and 25) where 
possible. The values are believed to be accurate to within 
one unit in the last place. 

For convenience of computation, the coefi&cients have been 
renumbered so that the pivotal points are coimted from the 
leading to the trailing edge. This renumbering is indicated 
by using indices in, s) instead of (m, p). The above results 
are formally imchanged, except that the roles of leading and 
trailing edge are interchanged in equation (C24), the sign of 
equation (C26) is reversed, and the sign of dp is reversed 
(eq. (C30)). The renumbered influence coefl&cients are 
given in tables I to TV. 
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TABLE I.— PIVOTAL POINTS, INCLINED PLAT PLATE SOLUTION, AND INFLUENCE COEFFICIENTS FOR p, X, o 


r» 

(for 

/" 

(for X) 

6.0273 

62.6483 

-2.4142 

-13.6669 

1.4066 

6.4797 

-1.0000 

-4.0000 

.6682 

2.8929 

-.4142 

-Z3431 

.1988 

2.0791 


Xn 

V Xm 

Tm 

(for 

a00961 

10.1632 

10.1632 

.03806 

6.0273 

-6.0273 

.08426 

3.2066 

3.2066 

.14646 

2.4142 

-2.4142 

.22222 

L8709 

1.8709 

.30886 

1.4966 

-1.4966 

.40246 

1.2186 

1.2185 

.60000 

LOOOO 

-LOOOO 

.69764 

.8207 

.8207 

.69134 

.6682 

-.6682 

.77778 

.6346 

.6345 

.86336 

.4142 

—.4142 

.01674 

.8033 

.8033 

.96194 

.1980 

-. 1989 

.99030 

.0985 

.0985 



TABLE n.— INFLUENCE COEFFICIENTS FOR VELOCITY (a) N=S 



Cm (velocity due to tblclmeiss) 


dM (velocity doe to camlDer) 



1 I 20.0060 -4.0719 0 -0.2242 0 -0.0719 0 

-7.6240 1L3137 -3.3692 0 -.2977 0 -.1329 

0 -4.3890 8.6601 -3.1543 0 -.3890 0 

-.6863 0 -3.4142 8.0000 -3.4142 0 -.6858 

0 -.3890 0 -3.1613 ^ 6691 —4.3890 0 

-.1329 0 -.2977 0 -3.3692 1L3137 -7.6240 


-.0719 0 -.2242 0 


-0719 20.9050 


0.2813 -0.0719 0.6791 

-.5700 0 -.8120 

.3914 -.2461 .0460 

-3. 6966 0 -1. 6307 

9.6368 -3.7549 2.1165 

-3.7606 11.3137 -9.6406 

7.1097 4.0719 38.0694 


(b) iV=16 


Cm (velocity doe to tblokness) 



2 

3 

-16.0614 

0 

41.8100 

-1L2032 

-16. 2646 

28,7992 

0 

-11.4300 

-L6324 

0 

0 

-1.1468 

-.4616 

0 

0 

3816 

-. 1958 

0 

0 

-. 1667 

-.0960 

0 

0 

-.0747 

-.0461 

0 

0 

-.0310 

-.0134 

0 



-1^6079 

41.8100 

-17.3778 

0 

-1.9240 

0 

-.4616 

0 

.1968 

0 

1.0183 

0 

3.6016 

0 

33.6568 




6 

B 

8 

9 

-0.1383 

0 

-0.0607 

0 

0 

-.1801 

0 

-.0764 

-.6896 

0 

-.2009 

0 

0 

-.6743 

0 

-.2172 

-6,9646 

0 

-.6735 

0 

17.3183 

-6.6633 

0 

6919 

-6.9676 

16.3135 

-6.4423 

0 

0 

-6.6685 

16.0000 

-6.6686 

-.7346 

0 

-6.4423 

16.3136 

0 

-.6919 

0 

-6.6633 

-.2361 

0 

—.6736 

0 

0 

2172 

0 

6743 

-.0942 

0 

-.2009 

0 

0 

-.0764 

0 

-.1801 

-.0262 

0 

-.0607 

0 


0.1518 
-L0788 
.1642 
-7. 9498 
10.4363 
-7.9207 
.2616 
-.9741 
.3736 
-.4142 
.6766 
-.2627 
L7841 
-.2079 
16.6480 


dm» (velocity doe to camber) 


-0.1337 

0 

-.2261 

0 

—.0492 

0 

-6.3694 

16.0000 

-6.2680 

0 

-.2474 

0 

1.1218 

0 

12.7610 






13 

13 

14 

15 

-0.0171 

0 


0 

0 

-.0310 

0 

-.0264 

-.0687 

0 

-.0461 

0 

0 

-.0747 

0 

-.0620 

-. 1804 

0 

-.0960 

0 

0 

-. 1667 

0 

-.1242 

-.3012 

0 

-. 1968 

0 

0 

-.3616 

0 

-.2699 

-.9353 

0 

-.4616 

0 

0 

-L1468 

0 

-.6457 

-00622 

0 

-L6324 

0 

22.6274 

-U.4300 

0 

-2.3602 

-8,9804 

28.7992 

-16.2646 

0 

0 

-11,2032 

41.8100 

-29.5439 

-.6512 

0 

-16.0614 

82.0133 


-0.0638 

0. 1618 

-00317 

0.2272 

1 

-0.0134 

06469 

0 

-.2079 

0 

-.2682 

0 

-.0733 

-.0942 

.1642 

-.0441 

.2467 

-.0182 

.6997 

0 

-.2627 

0 

-.3018 

0 

-.7617 

-. 1988 

.1033 

-.0816 

.2893 

-.0317 

.8238 

0 

-.4142 

0 

-.4024 

Q 

-.0480 

-.6447 

.2616 

-.1984 

.3766 

-.0692 

L0723 

0 

-.9741 

0 

-.6609 

0 

-1.3322 

-6.7791 

,3738 

-.7233 

.6500 

-. 1967 

1.5920 

17.3183 

-7.9207 

0 

-1.5233 

0 

-02377 

-6.6196 

10.0196 

-06626 

1.0125 

-.8740 

08834 

0 

-7. 9498 

22.6274 

-11.9890 

0 

-02438 

.6806 

1.7841 

-7.3218 

3L 4693 

-108582 

7.6037 

0 

-1.0788 

0 

-12.2147 

4L8100 

-37. 1476 

13.6103 

16.6480 

16.9635 

204100 

15.0614 

148.7048 
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EEPOET 1274 — ^NATIONAL ADVISOET COMMITTEE FOE AEEONAtmCS 


TABLE III.— INFLUENCE COEFFICIENTS FOR SLOPE 

(a) AT=8 


c— (slope doe to thickness) 


fnt (slope duo to camber) 


L5307 -0.8284 a6840 -0.6638 L0S24 

-4.7183 2.0000 -L4046 1.4142 -2.2068 


Z6131 -Z3978 3.6066 

-6.6568 4.0000 -6.6568 


-.4483 -4.8284 


2.3978 -2.6131 


L4046 -ZOOOO 


.4483 -8.0647 8.9218 

4.7183 L 4142 -17.0479 


.8284 -L6307 4.9932 


a3086 -9.2262 


3.6965 -2.1648 


L4142 -6.1048 


2.1648 -2.8284 


-L6307 1.8362 -2.6131 


1.2262 -1.4142 


1.6307 -1.2262 


2.8284 -1.8362 


.4483 -6.2262 


1.8362 -2.8284 


2.6131 -1.8362 1 

-6.2262 2.8284 -2 


-.4483 -0.1048 
6.1648 -1.4142 


1.2262 -1.6307 


2. 1018 -3. 0966 0. 2202 -0 


(b) iV^-16 
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TABLE IV.— INFLUENCE COEFFICIENTS FOR SECOND DERIVATIVE 

(a) N=8 


(b) JV=16 


;»• (second derivative doe to thlrimess) 


!■■■ 

2 

3 

4 

6 

6 


8 

9 

10 

11 

12 

13 

14 

16 

HI 
















1 

-6276.99 

1486.60 

-16L34 

36.17 

-1X16 

X36 

-X78 

1.63 

-101 

663 



123 

-X63 

5X47 

2 

1292.86 

-3102. 47 

72661 

-106 26 

3L 16 

-1X69 

631 

-X69 

X21 

-137 


.25 

-X44 

1131 

-101.83 

3 

044.06 

963.60 

-107X62 

44X16 

-78.01 

2X29 

-11.91 

X43 

-X86 

X88 

-123 

-.23 

X59 

-17.10 

15X04 

4 

-977,99 

-99.48 

648.07 

-67X00 

314.31 

-6X41 

2X14 

-1131 

X44 

-X90 

XIO 

0 

-161 

2X03 

-21117 

5 

844,68 


-9684 

366 76 

-489.31 

26X21 

-6110 

2165 

-1112 

646 

-X69 

.63 

X40 

-2X12 

2n.60 

6 

-710. 24 

4663 

2&13 

-7X63 

27X74 

-397.83 

31X 31 

-6X47 

2126 

-11.31 

X29 

-X02 

-X73 

3X31 

-337.63 

7 

604.40 


-4.08 

27.00 

-6X10 

23L18 

-36X29 

20X15 

-6X47 

22.15 

-1170 

194 

X13 

-4142 

41161 

8 

-496.08 

4688 

-X60 

-1L31 

04.21 

-5106 


-34X00 

21X07 

-6106 

24.21 

-1131 


4X88 

-49X08 

9 

411.61 

-4L42 

6.13 

4.94 

-11.70 

2X16 


20X15 

-35X29 

23118 

-62. 10 


-198 

-5X 14 

59140 

10 

-337.63 

3631 

-6.73 

-602 

X29 

-1L31 

2L26 

-60.47 

21X21 

-397.63 

27674 

-7X63 

2613 

4663 

-n62i 

11 

27L60 

-29.12 

640 

.63 

-669 

X46 

-U.12 

2L66 

-64.10 

266 21 

-4S9L21 

38678 

-96 84 

-2673 

844.68 

13 

-21L17 


-4.61 

0 

XIO 

-3.90 

X44 

-1131 

23.14 

-62.41 

814.31 


61607 

-99.48 

-977.99 

13 

166.01 

-17. 10 

669 

-.23 

-L23 

X38 

-686 

643 

-11.91 

2629 

-78.01 

442L16 


96660 

944.06 

H 

-101.83 

11.31 

-X44 

.25 

.68 

-L37 

X21 

-X69 

X31 

-12.69 

3116 


72X61 

-2162.47 

1202.86 

16 

60.47 


L23 

-. 16 

-.30 

.63 

-LOl 

163 

-X78 

636 

-12.16 

3X17 

-16134 

1486 60 

-627X 99 

Au (second derivative dne to camber) 

1 


2027.02 

-28X66 

8614 

-3669 

19.36 

-1L93 

X33 

-X48 

X62 

-X48 

616 

-638 

1646 

-64.63 

3 

4281. S3 

-2608. 18 


-144. 16 

6X26 

-24.00 

1X96 

-9.37 

7.11 

-X05 

X83 

-648 

X73 

^^snsjB 

6632 

3 


1206 89 


47L77 

-83.46 

3X72 

-1X48 

1167 

-638 

X90 

-X48 

7.06 

-X37 

1X97 

-69.31 

4 

483.46 

-28666 



337.17 

-7X48 

2X93 

-1600 


-X38 

7.66 

-600 

^^BnSsB 

-1648 

6X88 

5 

-286,66 

137.49 

-144.16 

885.71 

-60617 

267.28 

-69.31 

2X92 

-1X28 

1100 

-X37 

9.50 

-1193 


-7X48 

0 

194.70 

-73.94 


-03. 46 

28X34 

-40X83 

22X01 

-6163 

2X37 

-16 00 

IX 61 

-1193 



7X90 

7 

-144.16 

49.69 

-36.69 

4L66 

-7X48 

23X88 

-367.78 

210. 19 

-5163 

27.07 

-1X48 



2X94 

-9X45 

8 

113.41 

-30.09 



3L70 

-6X 31 

21136 

-34100 

21136 

-69.31 

3170 


2153 

-3X69 

11X41 

0 

-93.46 

28.94 

-18.09 

^^BrvuB 

-1X48 

27,07 

-6163 

21X19 

-357.78 

23X89 

-7X48 

4156 

-3X69 

4X69 

-14116 


7690 

-24.00 


-1L93 

1X61 

-1600 

2X37 

-5163 


-40X83 

28X24 

-9X46 

6X73 

-7X94 

19170 

11 


2a 71 



-9.37 

1100 

-1628 

2X93 

-69.81 

257.28 

-50X17 

38X71 

-14116 

127.49 

-28X66 

13 

03.88 

-ia48 



7.66 

-X38 

10.73 

-1X00 

2X93 

-7X48 

327.17 


60X90 

-28X66 

48X45 

13 

-60.31 

iao7 

-0.37 

7.06 

-X48 

X90 

-X38 

1167 

-1X48 

3X73 

-9X45 

47177 

-114X53 



14 

66.33 

-laoo 

8.73 

-6 48 

X83 

-X06 

7.11 

-9.37 

1X95 

-2100 

5X25 

-141 15 

82X31 


428183 

16 


15.46 

-&38 

a 16 

-X48 

X63 

-648 

X32 

-1193 

1X36 

-3X69 

8X14 

-28X06 


-1160153 


Am (second derivative dne to camber) 


1 

2 

3 

4 

5 

6 

7 

1 

-760.098 

144.162 

-21000 

9.373 

-6069 

X476 

-1X000 

3 

28X 658 

-181000 

70.481 

-16000 

X382 

-XOOO 

1X479 

3 

-73,041 

9X446 

-103.902 

61627 

-16000 

11927 

-21000 

4 

36686 

-101000 

59.314 

-88.000 

5X314 

-101000 

3X636 

5 

-21000 

11937 

-16000 

64.627 

-103.902 

93.446 

-7X941 

0 

IX 479 

-XOOO 

6382 

-1X000 

76481 

-181000 

28X 653 

7 

-1X000 

X476 

-6069 

9.373 

-21000 

141162 

-76X 008 


ffM (second derivative due to thickness) 


— 

1 

2 

3 

4 

5 

6 

7 

1 

-411392 

166260 

-12.686 

6687 

-1.373 

-0.260 

11.314 

a 

09.478 

-166000 

62.406 

-11314 

X806 

0 

-2XC32 

3 

46627 

7X529 

-97.608 

50.469 

-11314 

X019 

3X314 

4 

-46882 

-11.314 

K066 

-84.000 

51050 

-11314 

-4X882 

6 

36314 

2.019 

-11.314 

50.469 

-97.608 

7X629 

4X627 

6 

-23.032 

0 

X806 

-11.314 

62.406 

-16X000 

99.478 

7 

11314 

-.250 

-1373 

3.687 

-12.686 

106280 

-411392 


TABLE V.— COMPRESSIBILITY FACTORS (-y= 
















































































